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Abstract 



We introduce a fractional Kramers equation for a particle interacting with 
a thermal heat bath and external non-linear force field. For the force free case 
the velocity damping follows the Mittag-Leffler relaxation and the diffusion 
is enhanced. The equation obeys the generalized Einstein relation, and its 
stationary solution is the Boltzmann distribution. Our results are compared 
to previous results on enhanced Levy type of diffusion derived from stochastic 
collision models. 
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I. INTRODUCTION 



The pioneering work of Scher and Montroll [|IJ and Scher and Lax on the continuous 
time random walk || applied to diffusion problems led to a revolution in our understanding of 
anomalous diffusion processes. Anomalous diffusion is now a well established phenomenon, 
found in a broad range of fields 0-0]. It is characterized by a mean square displacement 

(x 2 )~t 5 (1.1) 

with (5^1. Various mechanisms are known to lead to enhanced diffusion 5 > 1 or sub- 
diffusion 5 < 1. Usually such processes are non-Gaussian meaning that the standard central 
limit theorem cannot be used to analyze the long time behavior of these phenomena. In order 
to describe such anomalous processes, fractional kinetic equations were recently introduced 
by several authors [|5|HTjJ- Within this approach, fractional space and/or time derivatives 
replace the ordinary time and/or space derivatives in the standard kinetic equation (e.g., 
Fokker-Planck equation). Examples include kinetics of viscoelastic media |§, Levy flights in 
random environments [10], chaotic Hamiltonian dynamics |Q|, and Quantum Levy processes 



|16| . For a discussion on Levy statistics and continuous time random walk in the context of 
single molecule spectroscopy see [|T7|-|T9||. 

In this paper, we introduce a fractional Kramers equation describing both the velocity v 
and coordinate x of a particle exhibiting anomalous diffusion in an external force field F(x). 
In the absence of the external force field the equation describes enhanced diffusion. The 
new equation we propose is 

dP{x,v,t) dP F{x)dP{x,v,t) 
dt dx M dv 

la0 Dl- a L fp P(x,v,t), (1.2) 

with < a < 1 and 

f 9 k b T d 2 

the dimensionless Fokker-Planck operator. Here we employ the fractional Liouville-Riemann 
operator [2(J oDl~ a in Eq. (|1.2| ) which we define later in Eq. (|2.4| ). j a is a damping 



coefficient whose units are [l/sec a ]. In the presence of a bounding force F(x) = —V'(x) 
the stationary solution of the fractional Kramers equation (|1.2j ) is the Maxwell-Boltzmann 
distribution. When a = 1 we recover the standard Kramers equation. 

Fractional kinetic equations are related to Montroll- Weiss continuous time random walk 
(CTRW) Il2~l~l-f2~4j1 . Here we show that the fractional Kramers equation is related to the 



coupled CTRW in the enhanced diffusion regime 5 > 1. This case corresponds to the so 
called Levy walks that are observed in a number of systems[4,5]. The different limits of the 
CTRW were used to model diverse physical phenomenon (when F(x) = 0) and therefore 
fractional kinetic equations in general and the fractional Kramers equation in particular are 
believed to be of physical significance. 

The basis of the fractional Kramers equation is the fractional Ornstein-Uhlenbeck process 
described by the fractional Fokker-Planck equation 
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?Q^=7aoDl-°itoQ(v,t), (1.4) 

Q(v, t) the probability density of finding the particle at time t with velocity v when F(x) = 0. 
We see that the fractional Kramers equation is a natural extension of Eq. (|1.4| ) in which the 
streaming terms describing Newtonian evolution are added in the standard way (i.e., as in the 
Boltzmann or Liouville equations). When a = 1 we get the standard Ornstein-Uhlenbeck 
process which has a fundamental role in non equilibrium statistical mechanics. 



The Rayleigh model was used to derive Eq. ( |1.4| ) for the normal case of a = 1 p5| , |26 
The Rayleigh model for Brownian motion, also called the Rayleigh piston, considers a one- 
dimensional heavy particle of mass M colliding with light non-interacting bath particles 
of mass m which are always thermalized. According to this model the moments of time 
intervals between collision events are finite. What happens when the variance of time inter- 
vals between collision events diverges? In this case we anticipate a non Gaussian behavior. 



This case has been investigated by Barkai and Fleurov |27| , g8| and as we shall show in 



subsection [11 B|, for a non-stationary model, such an anomalous case corresponds to the frac- 



tional Fokker-Planck equation (|1.4j) . However, as we shall show, the usual Rayleigh limit 
m/M — > in the non-Gaussian case is not as straightforward as for the ordinary Gaussian 
case. 

The more fundamental question of the necessary conditions for transport to be described 
by diverging variance of time intervals between collision events, or how to derive a from first 
principles is not addressed in this paper. In this context we note that several mechanisms 
in which the variance of time between collisions (or turning events in random walk schemes) 
diverge and which lead to anomalous type of diffusion are known in the literature 

Previous approaches |8|-15] have considered fractional kinetic equations in which the 



coordinate and/or time acquire the fractional character while in our approach the velocities 
are the variables that acquire fractional character. For F(x) = we find 1 < 5 = 2 — a < 2. 
The lower bound 5 = 1 corresponds to normal diffusion, the upper bound 5 = 2 can also be 
easily understood. For a system close to thermal equilibrium we expect 

(x 2 ) < (k b T/M)t 2 } (1.5) 

and hence 5 < 2. 

Previously Kusnezov, Bulgac and Do Dang (JH3] have suggested a fractional Kramers 
equation which was derived for the classical limit of a Quantum Levy process. Also in this 
work do the velocities acquire fractional character; however this approach is very different 
than ours, since it is based on Reisz fractional operators and gives (x 2 (t)) = oo (i.e., 5 = oo). 
Recently Metzler and Klafter considered a fractional kinetic equation that they also 
called a fractional Kramers equation. Their equation describes sub-diffusion (5 < 1) and is 
very different from our equation which described enhanced diffusion (5 > 1). 

This paper is organized as follows. In section [TJ we introduce the fractional Ornstein- 
Uhlenbeck process described by Eq. (|1.4| ), a brief introduction to fractional calculus is given. 



In subsection |11 A| the solution of the fractional Fokker-Planck Eq. (|1.4|) is presented and 



in subsection II B we derive this equation based on stochastic collision model. In section |Tl 1 
we consider the fractional Kramers equation ( |1.2| ), general properties of this equation are 
given and the force free case is investigated in some detail. We end the paper with a short 



summary in section IV 
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II. FRACTIONAL ORNSTEIN UHLENBECK PROCESS 



Let Q(v,t) be the probability density describing the velocity v of a macroscopic Brownian 
particle, with mass M, interacting with a thermal heat bath. The Fokker-Planck equation 
describing the time evolution of Q(v , t) with initial conditions Q(v, t — 0) is given by p6|j30 



d -^^ = ll L fp Q^t). (2.1) 

We shall call such a Fokker-Planck equation standard or ordinary. According to Eq. ( [2.1| ) 
the damping law for the averaged velocity is linear (v(t)) = —Ji(v(t)), and the velocity 
fluctuations are thermal. The stationary solution of Eq. (|2 . 1|) is the Maxwell's density 
defined with the thermal energy k^T. The corresponding Langevin equation is 

V = -Yiv + Z(t) (2.2) 



and £(£) is Gaussian white noise p6| , |30| . The stochastic process described by the Langevin 
equation is the well known Ornstein-Uhlenbeck process. Originally Eq. (|2.1| ) was derived by 
Rayleigh [^5j for a particle interacting with a gas consisting of light particles, however the 
scope of Eq. ( [2.1|) is much wider. It is used to model Brownian motion in dense environments 
when memory effects are negligible. 

In this section we generalize the Fokker-Planck equation ( |2.1| ) using fractional calculus. 
First we give some mathematical definitions and tools. 

The Liouville-Riemann fractional integral operator fl2(J of order a > is defined by 



* (t - t') a ~ ] 



D; a f(t) = / v ' f(t')dt'. (2.3) 
Jo 1 (a) 

For integer values a = n; oD^ n is the Riemann integral operator of order n. Fractional 
differentiation of order a > is defined by 

d n r i 

D?f(t) = — { Dr n f(t)}, (2.4) 

where n — 1 < a < n and d n /dt n is ordinary differentiation of order n. Within this fractional 
calculus 

Df^ = r , r( ^ + ^ n t^ (2.5) 
T (fji =f a + 1) 

when n > — 1. Notice that oDf a l ~ t Ta when < a < 1. The Laplace transform 

/•CO 

/(«)=/ e~ ut f(t)dt = C[f(t)] (2.6) 
Jo 

of a fractional Liouville-Riemann operator is 

n-l 

C Wf{t)\ = u a f(u) -Y,u k Dr l - k f(t)\t=o. (2.7) 



k=0 
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n is an integer satisfying n — 1 < a < n. For fractional integrals a < the sum on the RHS 
of Eq. (|2.7p vanishes. From Eq. (|2.7|) we see that the Laplace transform of a fractional 
derivative of f(t) depends on fractional derivatives of that function at time t = 0. In this 
work we use the convention that the arguments of a function indicate the space in which 
the function is defined, e.g. the Laplace transform of Q(v,t) is Q(v,u). 
The Fokker-Planck equation ([2.1|) is rewritten in the integral form 



Q(v, t) - Q(v, t = 0) = 7l Q Dt X Lf p Q{v, t) (2.8) 



We now replace the integer integral operator oD t 1 in Eq. ( |2.8j ) with a fractional integral 
operator oD^ a and < a < 1. The result is 

Q(v,t) - Q(v,t = 0) = la D- a L fp Q(v,t) (2.9) 

where j a is a generalized damping coefficient with units [sec] _a . Ordinary differentiation of 
Eq. ( |2.9| ) gives the fractional Fokker-Planck equation 



dQ(v,t) 
dt 



7a oD 



1-a 
t 



d k b T d 2 



dv M dv 2 



Q(v,t). (2.10) 



When a = 1 the ordinary Fokker-Planck Eq. ( [2.1|) is obtained. In Eq. fl2.10| ) we use the 
natural boundary conditions 

lim Q(v,t)= lim dQ(v,t)/dv = 0. (2.11) 

v— >±oo v— >±oo 

Later we shall show that the solution of Eq. ( |2.1U| ) is non-negative and normalized. Eq. 
( p,10| ) describes a fractional Ornstein-Uhlenbeck process. When the velocity v is replaced 
with a coordinate x, the equation describes an anomalously over damped harmonic oscillator 
as investigated in p3| . 

Eqs. (|2.9| ) and ( p.lOj ) are initial value problems. While Eq. (f2~9l) depends on a single 
initial condition [Q(v,t = 0) on the LHS of the equation], in solving Eq. (|2.10j ) two initial 
conditions have to be specified |20[| , these being Q(v,t = 0) and Q D^ a Q(v,t)\ t=Q . When 
setting D^~ a Q(v, t)\ t=Q equal to zero, the two equations are equivalent [51] . 

Multiplying Eq. ( |2.10| ) by v and integrating over v we find that the mean velocity is 
described by 

(v(t)) = - laQ D}- a (v(t)). (2.12) 
In Laplace space Eq. (|2.12|) reads 



U + 7aU l 

Vq is the initial velocity. The inverse Laplace transform of Eq. ([2.1 3|) is 

(v(t)) = v E a (- la r), (2.14) 

and 
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n=0 



r (1 + an) 



(2.15) 



is the Mittag-Leffler function. When a = 1 the Mittag-Leffler reduces to the exponential. 
For large t, Eq. Q2.14Q exhibits a power law decay 



(«(*)> 



^0 (7a*)" 



r(l-a) 

and for short times the relaxation is a stretched exponential 



(2.16) 



(v(t)) — v exp 
In a similar way we find the second moment 



T(l + a) 



(2.17) 



(v 2 (t)) = v 2 E a (-2 7 X) + -± f [l-E a (-2 7a t Q )] 



(2.18) 



exhibiting power law decay towards the thermal equilibrium value (v 2 (t = oo)) = k{I '/M. 
From Eqs. (|2.14 ) and (|2.18| ) we see that the Mittag-Leffler relaxation replaces the ordinary 
exponential relaxation found for ordinary Brownian motion. These equations were derived 
in pTl based upon a stochastic collision model which we will discuss in subsection |11B. 



A. Solution 



Our aims are (a) to find a solution of the fractional Fokker-Planck Eq. ( |2.10 ) and (b) 
show that Q(v, t) in Eq. (|2.10 ) is a probability density. The solution we find is an integral of 
a product of two well known functions. We use the initial conditions Q(v, t = 0) = 5{v — Vq) 
the generalization for other initial conditions is carried out in the usual way. 

We first show that the solution is normalized. The Laplace transform of Eq. (|2.10|) is 



1lQ(v, U) - S(V - V ) = "f a U a L fp Q(v,u). 



(2.19) 



Integrating Eq. ( 2.19 ) with respect to v, using the boundary conditions in Eq. (|2.1ip and 
the normalized initial condition we find 



f°° 1 
/ Q(v, u)dv = —. 

J-oo U 



[2.20) 



Since C(l,u) = 1/uwe see that Q(v,t) in Eq. ( [2.10|) is normalized. 
Let us now find the solution in Laplace u space. We write Q(v,u) as 



Q(v,u) 



where 



L fp G s (v) = —G 



R s (u)G s (v)ds 



d 



ds 



(2.21) 



(2.22) 
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and 

G (v) = 5(v-v ). (2.23) 
Eq. ( |2.22|) is the dimensionless ordinary Fokker-Planck Eq. ( |2.1| ), with solution p6| , |30 

a, M = 



M 



2irk h T (1-e 



-2s\ 



: exp 



M(v - v e~ s ) 2 
2k b T(l-e- 2s ) 



[2.24) 



We see that G s (v) is a non-negative probability density function, describing the standard 
Ornstein-Uhlenbeck process, and normalized according to 



G s (v)dv = 1. 



(2.25) 



R s (u) in Eq. ( |2.21| ) must satisfy a normalization condition. Using Eqs. ( |2.20|) and 
( p.25| ), we have 



oo I 

R s (u) ds = — . 

u 



(2.26) 



Inserting Eq. (|2.21| ) in Eq. (|2.19| ), using Eq. ( |2.22|) , and integrating by parts, we find 

u R s (u) G s (v)ds — 5(v — v ) = 



7aU 



.l-a 



,1-a 



R s (u) L fp G s (v) ds 

d 



R s (u) —G s (v) ds 

OS 



la u l ~ a [R^ (u) (v) - Ro (u) Go (v)} 



.l-a 



POO 


' d 




^Rs (u) 


lo 


OS 



G, (v ) ds. 



[2.27) 



According to Eq. (|2.26|) the boundary term ^^(m) in Eq. ( |2.27| ) is zero. Using the initial 
condition Eq. |E2|) in Eq. ( fL2l\ ) we find 



1°° \uRs (u)+ la u l - a 



d_ 

ds 



R. (u) 



G, (v) ds 



1 - 7 X~ a i2o(s) 5(v-v ). 



(2.28) 



Eq. ( |2.28| ) is solved once both sides of it are equal to zero; therefore, two conditions must 
be satisfied, the first being 



la u Rq (U) = 1 



(2.29) 
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and the second being 



7aU 



l-a 



d 

—R s (u) = uR s (u) . 

OS 



(2.30) 



The solution of Eq. fl2.30|) with the condition Eq. (|2.29| ) is 

1 



Rs (u) 



IolU 



l-a 



exp 



su 



To 



(2.31) 



It is easy to check that R s (u) is normalized according to Eq. ( 2.26|) . 

The solution of the problem in t space is the inverse Laplace of Eq. Q2.21 



Q{v,t) 



R s (t)G s (v)ds, 



where R s (t) is the inverse Laplace transform of R s (u) given by 

1 



R B (t) 



z a+1 L (z) 



(2.32) 



(2.33) 



and 



(la 



.I/O 



A/a 



(2.34) 



Properties of R s (t) are discussed by Saichev and Zaslavsky [|12| . l a (z) in Eq. (|2.33|) is one 
sided Levy stable density p2| , whose Laplace transform is 



l a (u) = / exp (—uz) l a (z) dz = exp(—u a ). 
Jo 



(2.35) 



The proof that R s (t) Eq. ([2.33| ) and R s (u) Eq.( |2.31|) are a Laplace pair is given in Appendix 
A (and see also 0]). 

A few features of the solution Eq. (|2.32j ) can now be discussed. An interpretation of 
Eq.( |2.32 ) in terms of a stochastic collision model will be given in the next subsection. 

When < a < 1; R s (t) is a probability density normalized according to / °° dsR s (t) = 1. 
Since G s (v) is also a probability density the solution Eq. ( |2.32|) is normalized and non 
negative. This justifies our interpretation of Q(v,t) as a probability density. 

When a = 1 the solution reduces to the well known solution of the ordinary Fokker- 
Planck equation. To see this, note that the inverse Laplace of Eq. (|2.31 ) for a = 1 is 
R s (t) = 1/jiS (t — s/jx), and then use definition Eq. ( |2.32| ). When a = 1/2 we have 



h/2(z) 



2A 



z 3//2 exp 



1 

Tz 



(2.36) 



with z > 0, then R s (t) Eq. ( [2.33|) is a one sided Gaussian 

Rs(t) = x __ 2 ^ exp 



\ nit/it 



H/2t . 



(2.37) 



S 



Two other closed forms of one sided Levy probability densities £2/3(2) and £1/3 (z) can be 
found in Series expansions of Levy stable density are in Feller's book |32] chapter 

XVII.6. 

In Fig. 1, we show the solution for the case a = 1/2 and for different times. The solution 
is found with numerical integration of ( |2.32| ) using Eqs. ( |2.24| ) and ( |2.37[ ). We choose the 



initial condition Q(v,t = 0) = 5(v — 1), (i.e., vq = 1) and k^T/M = 1. The solution exhibits 
a slow power law decay towards thermal equilibrium. We observe a cusp at v — Vq — 1; 
thus, initial conditions have a strong signature on the shape of Q(v,t). A close look at the 
figure shows that for short times (i.e., t < 2) the peak of Q(v,t) is at v — v Q . This is very 
different from Gaussian evolution for which the peak is always on (v(t)). 

Metzler, Barkai and Klafter |13] have shown that a fractional Fokker-Planck equation, 



which describes sub-diffusion 5 < 1, can be solved using an eigenfunction expansion which 



is identical to the ordinary expansion of the Fokker-Planck solution [|30|; but in which the 
exponential relaxation of eigenmodes is replaced with a Mittag-Leffler relaxation. We can 
use the eigenfunction expansion in |13| to find a second representation of Q(v , t) in terms of a 



sum of Hermite polynomials. Expanding G s (t>), using the standard eigenfunction technique 
of Fokker-Planck solutions, we write 



I ^ G/5*) (&) exp { - ns) (2 - 38) 



where H n are Hermite polynomials. We insert the expansion Eq. (|2.38|) in Eq. ( |2.32|) and 



use 



POO 

/ R s (t) exp (-ns) ds = E a (-n^ a t a ) , (2.39) 
Jo 

to find the eigen function expansion 



^. . M ( Mv z 

Q {v,t) = \ — — — exp 



2nk h T \ 2k h T 




n=o z n - 



-n la t a ) . (2.40) 



The stationary solution, determined by the smallest eigen value n = 0, is the Maxwell 
distribution which is independent of 7 a and a. 

An extension of the fractional Fokker-Planck equation ( |2.10| ) to higher dimensions is 
carried out by replacing the one- dimensional Fokker-Planck operator, Eq. (|1.2|) with the 
appropriate d dimensional Fokker-Planck operator (e.g., replace d/dv with V). The solution 
for such a d dimension equation is then found to be Eq. (|2.32|) in which Eq. ( 2.24 ) must 



be replaced with the appropriate solution of the d dimensional ordinary Fokker-Planck 
equation. 
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B. Collision Model 



As mentioned above, the ordinary Fokker-Planck equation Q2.1| ) was derived by Rayleigh 
over a century ago. Briefly, the Rayleigh model for Brownian motion considers a one dimen- 
sional test particle with mass M colliding with bath particles of mass m and e = m/M « 1. 



When collisions are frequent but weak, the ordinary Fokker-Planck equation is valid j26 
Here we consider the case when the concept of collision rate does not hold and the mean time 
intervals between collision events diverges. As mentioned in the introduction this case was 
investigated in [^7],^J and as we shall show now such a case corresponds to the fractional 



Fokker-Planck equation ( |2.1(J| ). 

We consider a particle of mass M which moves freely in one dimension and at random 
times it collides elastically with bath particles of mass m. Bath particles are assumed to be 
much faster than the test particle. Collisions are elastic and one-dimensional and therefore 
the velocity Vj^ of the test particle immediately after a collision event can be related to the 
velocity of the test particle just before the collision event according to 

= (^) Vm + ^f m . (2.41) 

where e = m/M and v m is the velocity of bath particle distributed according to Maxwell's 
distribution. 

The times between collision events are assumed to be independent identically distributed 
random variables implying that the number of collisions in a time interval (0, t) is a renewal 
process. This is reasonable when the bath particles thermalize very quickly and when the 
test particle is slow. According to these assumptions the times between collision events {r^} 
are described by a probability density if) (r) which is independent of the mechanical state 
of the test particle. Therefore the process is characterized by free motion with constant 
velocity for time T\ then a collision event described by Eq. ( [2.41 ) and then a free evolution 



for a period T2, then again a collision etc. The most important ingredient of the model is 
the assumption that if>{r) decays like a power law for long time 

V>(t) ~ r- x - a (2.42) 



with < a < 1. From Eq. ( 2.42| ) we learn that the mean time between collision diverges, 



/ °° rip(r)dr = oo. Since in such a problem there is no characteristic time scale, the number of 
collisions in an interval (0, t) is not proportional to t for large times. In other words the law of 
large numbers is not valid for the choice Eq. ([2 .42 ), leading to non normal behavior. Similar 



waiting times distributions were used within the CTRW to model anomalous diffusion for 
the past three decades. When if>(r) is exponential and in the presence of an external force 
field F(x), this model was investigated extensively in the context of reaction rate theory 

IH3- 

Such a model is non stationary and the probability of a collision event in a small time 
interval (t, t + dt) is time dependent even in the limit of large times. This is a consequence 
of the diverging first moment. 

Let Qcoi (v, t) be the probability density for finding the test particle with velocity v at 
time t and initially v — Vq. Using the model assumptions 
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Qcoi (v,t) = J^R s (t) G s (7 



(2.43) 



s=0 



with R s (t) is the probability that s collision events have occurred in the interval (0,i) and 
G s (v) is the conditional probability density of finding the particle with velocity v after s 
collision events. We note that Eq. ( |2.43| ) is the discrete version of Eq. ( |2.32| ) . 
Using the map Eq. ( |2.41| ), it can be shown that 



a (v) 



M 



2nk b T (1 - n\ 



2s 



: exp 



M (y - v fip 2 
2k b T(l-fij s ) 



(2.44) 



with fi = (1 — e)/(l + e). Not surprisingly G s {v) is Gaussian since velocities of colliding 
particles are Gaussian random variables. We note that 



with the initial condition 



dG s (v) f -i\ f r, ( n 
Qs = ln [Hi ) L fpGs [V ) 



G (v) = S (v - v ) 



(2.45) 



(2.46) 



and ln (/^r 1 ) ~ 2e. Eq. (|2.45|) is a Fokker-Planck equation in which s, the collision number, 
plays the role of dimensionless time. 

The Laplace t — > u transform of R s (t) , R s (u) can be calculated using renewal theory 



R* (u) 



1 - lj}(u) 



u 



and ip(u) = C[if)(r)]. From Eq. ( |2.43|) we have in Laplace t — > u space 

oo 

Qcol (V, U) = J2 R s ( U ) G s (v) , 



(2.47) 



(2.48) 



s=0 



multiplying this equation from the left with Lf p , using Eq.( |2.45| ), and integrating by parts, 
we have 

LfpQcoi {v, u) = — 2^ R s (u) 



s=0 



ds 



1 ~ 1 00 

Rq (u) 8 (V — Vn) ^ 

2e 0V ' V 07 2e^ 



G s (v) 



(2.49) 



We have used Eq. ( |2.46| ) and the the boundary condition R^ (u) = for u ^ 0. 
According to Eq. ( p. 47] ) 

^R s (u) = ln R s (u) 



(2.50) 



and since according to (|2.42|) 



r ~ r 



-l-a 



we have 
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^){u)~\-Au a (2.51) 

(here A is a parameter with units t a ) valid for small u, inserting Eq. (|2.51| ) in Eq. (p. 50|) 
we have 

d ~ 

—R s (u) Au a R s (u) (2.52) 

OS 

and from Eq. ( p. 471 ) 

R (u) ~ Au a ~\ (2.53) 

We are now ready to derive the fractional Fokker-Planck equation. Inserting Eqs. (|2.52|) 
and ( |2.53| ) in Eq. ( |2.49| ) we find in the limit of small Au a and e 

uQcoi (v, u)-5(v- v ) = "faU^LfpQcoi (v, u) (2.54) 

with 

7q = lim ^. (2.55) 

Eq. ( |2.54| ) is the fractional Fokker-Planck equation ( j2.1| ) in Laplace space. We note that 
the moments of the collision model in the limit e — * found in |27j are identical to those 
found here based upon the fractional Fokker-Planck equation Eqs. ( |2.14| ) and ( |2.18| ) as they 
should be. 



III. FRACTIONAL KRAMERS EQUATION 

Let P(x,v,t) be the joint probability density function describing both the position x 
and the velocity v of a Brownian particle subjected to an external force field F(x). The one 
dimensional Kramers equation models such stochastic motion according to 

dP(x,v,t) dP F{x)dP(x,v,t) 
dt dx M dv 

7iL /p P(s,w,t), (3.1) 

where the Fokker-Planck operator L f p is given in Eq. ( |1.3| ) . The Kramers Eq. ( |3.1| ) implies 
that noise is white and Gaussian and it describes under-damped motion close to thermal 
equilibrium. Eq. ( |3.1| ) is an extension of Eq. ( |2.1| ) which includes the coordinate x as well as 
the effects of F(x). We generalize Kramers equation in the same way as above, and consider 

dP(x, v, t) +v dP(x, v, t) F(x) dP(x, v, t) 
dt dx M dv 

la0 Dl~ a L fp P(x,v,t), (3.2) 

with < a < 1. The terms on the LHS of the equation are the standard streaming 
terms describing reversible dynamics according to Newton's second law of motion. The 
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term on the RHS of the equation describes an interaction with a bath, it can be considered 
as a generalized collision operator replacing the ordinary collision operator found in the 
standard Fokker-Planck equation. As mentioned in the introduction the stationary solution 
of Eq. ( |3.2|) is the Maxwell-Boltzmann distribution and when a = 1 we recover the ordinary 
Kramers equation. 

A formal solution of the fractional Kramers equation can be found in terms of the solution 
of the ordinary Kramers equation. We denote the solution of the fractional Kramers equation 
with P a (sc, v, t, 7 a ) instead of P(x,v,t) we have used so far. The Laplace transform of Eq. 
(Qis 

uP a (x, v, u, ja) - P a (x, v,t = 0, 7„) 

| ^ dPg (x,v,u,jq) F(x) dP a (x,v,u,j a ) _ 
dx M dv 

n/ a u 1 ~ a L fp P a (x,v,u,7 a ), (3.3) 



and P a (x, v , u, j a ) is the Laplace transform of P a (x, v, t, j a )- From Eq. (|3.3| ) we learn that 
P a {x, v, u, 7 a ) solves an ordinary Kramers equation in which the damping coefficient 71 was 
transformed according to 

7i -> lau 1 -*. (3.4) 

We therefore find 

P a (x, v, u, ja) = Pi (x, v, u, 7 a M 1_ °) , (3.5) 

assuming the initial conditions are identical for both solutions. Transforming to the time 
domain we find the formal solution 

P a (x, v,t,j a ) = C' 1 Pi (x,f,M,7 Q M 1_Q ) (3.6) 

with £ -1 being the inverse Laplace transform. Closed form solutions of ordinary Kramers 
equation, P\ (x, v, t, 71), are known only for a handful of cases, approximate solutions can be 
found using methods specified in [[30[| . In some cases Eq. ( |3.6|) can be used to find moments 



of the solution of fractional Kramers equation, (x n v m } } in a straight forward way. In what 
follows we shall revert to the notation P(x, v, t) instead of P a (x, v, t, 7 a ). 



A. Force free case 1 

We consider the force free case F(x) = 0. As usual (x(t)) = (v(t)), with the mean 
velocity (v(t)) given in Eq. ( |2.14j) , hence 

(x(t)) = v tE a>2 (- la t a ) (3.7) 

where 

00 z k 
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is a generalized Mittag-Leffler function [|!| satisfying 



N-l 



with z — > oo. For short times 



E 



j r (/3 — an) 



+ O (z~ n ) 



(x(t)} ~ u t 

as expected for a pure ballistic propagation, and for long times 



(x(t)) 



v t 



1-a 



laT (2 - a) 



(3.9) 



(3.10) 



(3.11) 



The particle exhibits a net drift in a direction determined by the initial velocity vq. Of 
course, when averaging over initial conditions, using thermal equilibrium condition, no net 
drift is observed as expected from symmetry. The mean square displacement is determined 
by (x 2 (t)) = (2x(t)v(t)). A short calculation using the Laplace transform of Eq. ( |3.2|) shows 



<* 2 (t)> 



eq 



(3.12) 



where the subscript eq means that thermal initial conditions are considered (i.e., (v 2 ) 
kbT/M). For short times 



(x\t)) 



eq 



foT 2 

~ 1 , 



M 



while for long times 



where 



(x\t)) 



eq ~ 2D a t 



2~a 



k h T 



7«Mr (3 - a) ' 



(3.13) 



(3.14) 



(3.15) 



Eq. (|3.14j ) exhibits an enhanced diffusion when < a < 1. Eq. (|3.15| ) is the (first) 
generalized Einstein relation and when a = 1 we recover the well known Einstein relation 
D 1 = (fc 6 T)/(M 7l ). 

It is straightforward to prove the more general Einstein relation 



(x 2 (t)) 



eq 



dt' / (v(r)v(0)) eq dr, 



where according to Eq. (|2.14|) 



(v(t)v(O)) 



eq 



M 



(3.16) 



(3.17) 



is the velocity autocorrelation function. We note that relation Eq. ( |3.16j ) is valid for sta- 
tionary processes [|3|, while the collision model we investigated in previous section is non- 
stationary. 
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Now consider the constant force field F(x) = F, using the fractional Kramers Eq. fl3.2|) . 
We can show that a second generalized Einstein relation between the drift in the presence 
of the driving field, (x(t))p, and the mean square displacement, Eq. ( |3.12| ), in the absence 
of the field is valid 

«t)) F = F { ^^. (3.18) 

This relation suggests that the fractional Kramers equation is compatible with linear re- 
sponse theory ■ 



B. Force free case 2 



We shall now use the formal solution Eq. (|3.6|) to find moments (x 2n (t)) eq for the case 
F(x) = 0. Odd moments are equal zero. Consider the reduced probability density W eq (x,t) 
of finding the particle on x at time t, defined according to 



W eq {x,t) 



dv / dv P(x,v 1 t)M(v ) 



(3.19) 



P(x,v,t) is the solution of the fractional Kramers equation with initial conditions concen- 
trated on xq and vq. M(vq) is Maxwell's probability density implying an equilibrium initial 
condition for the initial velocity Vq. 

For the standard case a — 1, W eq (x, t) is a Gaussian, therefore 



(x 



2n 



(*)> 



(:'(./ 



x W eq (x, t) dx 



(2n)! 
2 n n! 



eq 



(3.20) 



and according to ^0 



(x 2 (t)> 



k b T[yit- l + exp(-7it)] 



eq 



M 



7? 



which in the long time limit gives (only for a 



(3.21) 



(x 2 (t)) 



eq 



k b T 
My! 



(3.22) 



According to Eqs. (p.4| - ^6|) the calculation of (x 2n (t)) eq , for < a < 1, follows three 



steps. First find the Laplace transform of the Gaussian moments, (x 2n (t)) eq Eq. (3.20). 



Since we shall be interested in the long time behavior of (x 2n (t)) eq it is sufficient to consider 



only the long time behavior of the standard a = 1 case, namely we use the asymptotic Eq. 
( |3.22j ) instead of Eq. ( 3.21|) . We note that the inclusion of the short time behavior is also 
straightforward, but of less interest to us here. It is easy to show that 



(x 2n {u)) eq ~(2n)\ 



My x 



u 



n+l 



(3.23) 
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valid for small u and a = 1. The second step is to transform Eq. (|3.23|) using Eq. ( |3.4| ) and 
the last step is to invert Laplace transform the result from the second step, as in Eq. fl3.6|) . 
We find 



(x*(t)> e9 ~(2n)! 



k h T 



j.n(2-a) 



M>y a l r (2n - na + 1) ' 



(3.24) 



It is easy to check that the moments in ( p.24|) can also be calculated based on 



(x 2n (t)> ~ /r 1 



'A' 



2ra 



U 



l-a 



U 



2-a i k b T 



M 7q 



A- 2 



fc=0 



(3.25) 



This result is also to be expected. In Fourier-Laplace space the well known a = 1 solution 
is 



W eq {k, u) 



u+^k 2 ' 



(3.26) 



if we use the transformation ( |3.4|) in Eq. ( |3.26|) we find 

,1-Q 



W eq (k, u) 



u 



u 



2-a 



M7, 



(3.27) 



And this is the moment generating function in Eq. ( |3.25| ). We note that in general there is 
no guarantee that the transformation Eq. ( p.4|) can be made after the small u limit of the 
a = 1 solution is taken, instead the small u limit must be taken only after the transformation 
Eq. (|3.4| ) takes place. The small (k, u) expansion a particular coupled CTRW also known 
as Levy walks has the same form as Eq. ( |3.27| ) {see Eq. 38 in jR)]]} . 

According to Eq. ( p.27| ) W eq (x,t) satisfies the following fractional diffusion equation 



dW eq (x,t) _ k b T 



dt 



dx 2 



W eq (x,t) 



(3.28) 



which is expected to work well only at large times. Eq. ( |3.28|) was investigated by Schneider 
and Wyss 0. 

For one dimension the inverse Laplace-Fourier transform of Eq. (|3.27|) is 



W eq (x,t) 



VlaM 



V2^T(1 -a/2) \t 



T^afi ll ~ a / 2 W 



(3.29) 



with 



z = t 



Vk^r 



\x\ 



l-a/2 



(3.30) 



The inversion of the two and three dimensional solutions is not as straightforward as for the 
one dimensional case, and we leave the details for a future publication. 



16 



C. Collision Model 



We have derived the fractional Fokker-Planck equation (|2.1Up from a non- stationary 
stochastic collision model in the limit of e = m/M — > 0. Eq. ( [2.10| ) describing a fractional 



Ornstein-Uhlenbeck process, is the basis of the fractional Kramers equation which is reached 



after adding the standard streaming terms (i.e. Newton's law of motion). In [27] , 28| the mean 
square displacement, {x 2 (t)) co i of the test particle was calculated also for mass ratios e not 
tending to zero. It was shown the collision process gives (x 2 ) co i ~ t 2 for < a < 1 and finite 



e, a behavior different from the asymptotic behavior we have found in Eq. (|3.14j) based 



on the fractional Kramers equation, (x 2 ) ~ t 2 ~ a . The (x 2 ) co i ~ t 2 behavior of the collision 
model can be understood by the fact that within the collision process, one typically finds 
long (collision-free) time intervals, of the order of the observation time t, in which the motion 
is ballistic so that x 2 ~ t 2 . A similar ballistic behavior is known from the Levy walk model 
with diverging time intervals between turning points. As e becomes smaller, the correction 
to the ballistic term becomes important for longer times and these correction terms behave 
like (x 2 ) ~ t 2 ~ a . Also it is clear that the two limits of t — > oo and e — ► do not commute. 
If we first take t — > oo and only then e->0we find 

(x 2 (t)) col ~ t 2 . (3.31) 

Hence the derivation of the fractional Kramers equation from the stochastic collision model 
is a delicate matter, with a result depending on the order in which limits are taken. Another 
difficulty in our derivation is that our starting point is a non-stationary process. It is 
still unclear if stationary models can lead to dynamics described by the fractional Kramers 
equation. 



IV. SUMMARY AND DISCUSSION 



We have investigated a fractional Kramers equation which has the following properties: 

(a) the velocity of the particle evolves according to a fractional Ornstein-Uhlenbeck process 
described by Eq. ( |2.10| ), the velocity moments decay according to a Mittag-Lefner relax- 
ation, namely as a stretched exponential (Kohlrausch form) for short times and as a power 
law for long times, 

(b) in the absence of a force field diffusion is enhanced, 1 < 5 < 2, 

(c) the stationary solution of the fractional Kramers equation is the Maxwell-Boltzmann 
distribution, 

(d) Einstein relations are obeyed in consistency with fluctuation dissipation theorem and 
linear response theory, 

(e) in Laplace space a simple transformation of solutions of ordinary Kramers equation gives 
the solution of fractional Kramers equation. 



As mentioned in the introduction fractional kinetic equations in the literature are related 
to the CTRW. We showed here that in the small (k, u) limit W eq (k, u) has the same form as a 
particular coupled Levy walk CTRW. Other limits of the CTRW are shown to correspond to 
other fractional kinetic equations. As pointed out in [2^] the fractional diffusion equation || , 
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describing the sub-diffusion, corresponds to the uncoupled CTRW in the limit of small (k, u). 
A comparison between the uncoupled CTRW and solution of fractional diffusion equation 
in (x,t) space was carried out in ||24j| . Another fractional equation in fllUyilf describes Levy 
flights for which (x 2 ) = oo, such a fractional equation is related to the decoupled limit of 
the CTRW with diverging jump lengths. 

One may ask if it is worthwhile to introduce fractional derivatives, given that the older 
CTRW approach is so successful. Besides the fact that fractional equations are beautiful 
and simple (i.e., in some cases they are solvable), these equations can incorporate the effect 
of an external potential field. To us, this extension seems important although not explored 
in depth in the present paper. Little is known on anomalous diffusion in an external force 
field. 



Metzler, Barkai and Klafter [13j have investigated a fractional Fokker-Planck equation 
defined with the fractional Liouville -Riemann operator. In the absence of an external 
force the fractional Fokker-Planck equation investigated in |T3j describes a sub-diffusive 
behavior (5 < 1). The equation considers a type of over damped dynamics in which only the 
coordinate x is considered not the velocity v. The fractional Fokker-Planck equation in fi"3 



together with the fractional Kramers equation investigated here give a stochastic framework 
for both sub and enhanced diffusion in an external field. We believe that both approaches 
will find their application. 

Note added in proof. Recently related work on fractional diffusion was published j|7 . 
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VI. APPENDIX A 

We find the Laplace transform of Eq. ( |2.33| ) 



R,(u) 



-Lit 



-z a+1 l a (z) dt, 



(6.1) 



with z defined in Eq. ( |2.34| ). Using the change of variables t = y{s/ r y a ) 1 ' a it is easy to show 
(s / 0) 



1 d f°° 

R s (u) = — / exp 

as du Jo 



1 \ d 




- — — exp u 

,asj du \ 7 a 



L (y) dy 



u a 1 ( su c 
exp 



(6.2) 



which is R s (u), Eq. ( 12.311) . From Eq. (|2.29|) we learn that Eq. (|2.31| ) is valid also for s = 0. 
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Figure Caption 

Figure 1: The dynamics of Q(v, t) for the fractional Ornstein-Uhlenbeck process with a = 
1/2 and for times, t = 0.02,0.2,2, 20 (solid, dashed, dotted, and dot-dash lines, respectively). 
Also shown (fine dotted curve) is the stationary solution which is Maxwell's distribution. 
Notice the cusp on v = vq = 1 as well as the non-symmetrical shape of Q(v,t). 
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